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Abstract 

Causal  FIR  matrices  with  anticausal  FIR  inverses  have 
a  key  role  in  the  theory  of  FIR  perfect  reconstruction 
filter  banks.  We  explore  the  theory  of  such  matri¬ 
ces.  Some  general  results  on  nature  of  inverses  of  first 
order  caus^  FIR  matrices  are  then  presented.  This 
leads,  in  particular,  to  a  complete  parameterization 
of  the  Biorthonormal  Lapped  Transform  (BOLT)  re¬ 
ported  earlier  in  (1].  ^ 


The  new  analysis  and  synthesis  filters  are  z~‘**Hk(z) 
and  z^*^Fkiz),  which  means  that  the  frequency  char¬ 
acteristics  are  unchanged.  By  choosing  /  large  enough 
we  can  make  Eifz)  causal  and  its  inverse  Ri(z)  anti- 
causrd.  So  the  class  of  all  FIR  PR  filter  banls  can  be 
characterized  by  considering  the  cafacafi  class.  Unlike 
arbitrary  FIR  systems  with  FIR  inverses,  cafacafi  sys¬ 
tems  have  special  system-theoretic  properties  [2], [3], 
and  generally  speaking,  are  more  interesting.  In  [3j 
it  is  shown  that  the  most  general  degree  one  cafacafi 
system  has  the  form  GoV^fz)  where 


1.  INTRODUCTION 

It  is  well-known  that  the  Af-channel  maximally  dec¬ 
imated  analysis/synthesis  system  of  Fig.  1.1(a)  can 
be  redrawn  as  in  Fig.  1.1(b)  where  E(z)  and  R(2) 
are  the  polyphase  matrices  of  the  analysis  and  synthe¬ 
sis  bank  respectively.  In  a  recent  paper  [1]  we  argued 
that  the  set  of  all  FIR  perfect  reconstruction  (PR)  sys¬ 
tems  can  be  essentially  characterized  by  characterizing 
all  CAusal  Fir  transfer  matrices  with  AntiCAnsal  Fir 
Inverses,  abbreviated  cafacafi.  In  particular,  this  gave 
rise  to  the  so-called  BiOrthnormal  Lapped  Transform, 
abbreviated  BOLT.  Some  theorems  pertmning  to  the 
factorization  of  cafacafi  systems  and  BOLT  systems, 
as  well  as  the  most  general  degree  one  cafacafi  build¬ 
ing  block  were  presented  in  [1],  along  with  examples  of 
unfactorizable  cafacafi  systems.  After  a  brief  review  of 
the  results  in  |1],  we  present  new  results  for  first  order 
polyphase  matrices,  and  establish  a  complete  factor¬ 
ization  of  the  BOLT  (which  was  missing  in  [1]). 

1.1.  Importance  of  cafacafi  systems 

The  filter  bank  in  Fig.  1.1  is  a  perfect  reconstruc¬ 
tion  (PR)  system  [i.e.,  x(n)  =  i(n)]  if  and  only  if 
R(2)  =  E~^(2).  It  is  an  FIR  filter  bank  if  both  of 
the  matrices  E(z)  and  R(z)  are  FIR.  So,  for  an  FIR 
PR  system  we  should  have  [det  E(z)]  =  cz~'^,  c  0. 
In  fact  all  FIR  PR  systems  are  characterized  by  this 
property.  Given  an  arbitrary  FIR  PR  system,  if  we 
define  Ei(z)  =  2“^E(z)  and  Ri(z)  =  z^R(2)  for  arbi¬ 
trary  I,  we  have  constructed  a  new  FIR  PR  filter  bank. 
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V„(z)  =  I-UmvjL-bZ”^Umv]L,  (1.1) 

Go  is  nonsingular,  and  and  v,„  are  M  x  1  vec¬ 
tors  with  UmVm  =  1.  (Superscript  f  stands  for  trans¬ 
pose  conjugate).  With  =  Vm,  V,„(z)  is  the  famil¬ 
iar  degree-one  paraunitary  (PU)  buildi^  block  [4], [5]. 
It  is  clear  that  Vm(l)  =  I;  it  can  be  shown  that 
(det  Vm(2)]  =  z“*.  Inspite  of  the  generality  of  this 
form,  an  arbitrary  cafacafi  system  may  not  be  fac- 
torizable  into  these  building  blocks!  An  example  is 

V(z)  =  uvt  -I-  z~*I  +  z“*uv^,uW  =  0,  see  [3]. 

1.2.  The  BOLT 

A  BOLT  is  a  first  order  causal  FIR  system 

G(z)  =  g(0)  +  2~*g(l),  (1.2) 

with  anticausal  FIR  inverse.  In  short,  a  BOLT  is  a 
M  xM  first  order  cafacafi  s^tem.  Its  degree  is  clearly 
equal  to  the  rank  of  g(l)  [5]. 

Recall  that  the  lapped  orthogonal  transform  (LOT) 
[6],  (7],  is  essentially  a  filter  bank  in  which  the  polyphase 
matrix  E(z)  is  paraunitary  with  order  one.  This 
means  that  £(z)  has  the  form  (1.2),  and  its  inverse  is 

gt  (0)  -b  zg^  (1)  which  is  FIR  anticausal!  So  the  BOLT 
is  a  generalization  of  the  LOT;  the  LOT  is  paraunitary, 
whereas  the  BOLT  merely  cafacafi. 

It  turns  out  that,  while  arbitrary  cafacafi  systems 
cannot  be  factorized,  any  BOLT  can  be  factorizi^  into 
degree  one  building  blo^: 


G(z)  =  G(l)V^(z)V^_,(z) . . .  Vi(z)  (1.3) 


where  axe  as  in  (1.1)  with  UmV^  =  1  (Fig.  1.2). 

This  therefore  is  a  generalization  of  the  paraunitary 
factorization  (4)  for  the  first  order  cafacafi  case.  The 
PU  factorization,  however  holds  for  any  order. 

Given  that  each  building  block  Vm(zl  is  cafacafi, 
it  is  easily  verified  that  the  product  (1.3)  is  cafacafi. 
However,  the  product  in  general  is  not  a  BOLT,  that 
is  it  may  not  be  a  first  order  system  as  in  (1.2).  This, 
however,  can  be  ensured  by  imposing  further  condi¬ 
tions  on  Ui  and  Vj.  This  result  (Sec.  3)  will  follow 
directly  from  the  more  general  results  on  first  order 
systems  to  be  developed  in  the  next  section. 


2.  INVERSES  OF  FIRST  ORDER  SYSTEMS 


Consider  a  first  order  M  x  M  transfer  matrbc  as  in 
(1.2).  We  will  assume  that  G(l)  =  G(e^°)  is  nonsin¬ 
gular.  This  is  the  case,  for  example,  when  there  exists 
an  FIR  inverse.  We  can  then  express  G(z)  =  G(l)F(r) 
where  F(l)  =  I.  So  we  can  write  F(z)  =  I  —  P  -H  z~^P 
without  loss  of  generality.  Let  p  be  the  rank  of  P  [i.e., 
p  =  degree  of  F(2:))  so  that  we  can  write 

F(2)  =  I  -  -t- (2.1) 

where  Li  and  V  are  Mxp.  The  type  of  inverse  that  F(2) 

has  depends  strongly  on  the  matrix  V^U,  in  particular 
on  its  eigenvalues.  We  now  present  some  precise  results 
on  this.  Throughout  this  section  F(2)  is  the  M  x  M 
system  shown  above;  U  and  V  are  Mxp  with  rank  p. 

Lemma  2.1.  F(2)  has  an  anticausal  inverse  if  and 
only  if  V^U  is  nonsingular.  0 

Proof.  A  minimal  implementation  of  Ffz)  is  shown 
in  Fig.  2.1.  The  state  space  description  of  this  struc¬ 
ture  (with  standard  notations  [5])  is  given  by  A  =  0, 

B  =  Vt,  C  =  W,  and  D  =  I  -  Wvt  Let 


A  B 

0  vt  ■ 

C  D 

U  l-wvt. 

As  shown  in  [2],  7^  is  nonsingular  if  and  only  if  there 
exists  an  anticausal  inverse  for  F(2).  We  will  show  that 

TL  is  nonsingular  if  and  only  if  V^U  is  nonsingular. 
Suppose  Tix.  =  0  for  some  vector  x.  Then 

V^X2  =  0,  and  Wxi  -I-  X2  =  0,  with  x  = 


S 


Combining  these  we  get  V^^Wxi  =  0.  If  is  nonsin¬ 
gular  then  xi  =  0  and  so  X2  =  —Uxi  =  0  from  above. 
In  short,  if  Tlx.  =  0  then  x  =  0.  So  7^  is  nonsingu¬ 
lar.  On  the  other  hand,  if  V^U  is  singular  there  exists 
y  ^  0  such  that  V^Wy  =  0.  If  we  now  choose  xj  =  — y 
and  X2  =  Uy,  then  Tix  =  0  so  7i  singular.  VV  V 

Theorem  2.1.  The  inverse  of  F(2)  is 


2.  FIR  and  anticausal  (i.e.,  F(z)  is  cafacafi)  if  and 

only  if  UV^  has  p  of  its  eigenvalues  equal  to  unity 
and  the  remaining  M-p  eigenvalues  equal  to  zero. 

3.  FIR  and  causal  (i.e.,  F(z)  is  unimodular  in  z~^) 

if  and  only  if  UV^  has  all  eigenvalues  =  0.  0 

Proof.  FVom  the  unitary  triangularization  theo¬ 
rem  (8)  we  can  write  UV^  =  TAT^  where  TTt  = 
I,  and  A  is  upper  triangular  with  the  ei^nvaules 
{Ao,  Ai . . .  Ap_i, 0, ... 0}  on  the  diagonals.  (Since  the 
rank  is  p  there  could  be  at  most  p  nonzero  eigenvalues). 

We  can  then  express  F(z)  =  T(l  —  A  2~*A)Tt  so 
that  (det  F(z)]  =  nf=o  (l  ~  This  has  the 

form  cz~^  ^hich  is  necessary  and  sufficient  for  the 
existence  of  FIR  inverse)  if  and  only  if  Aj  =  0  or  1  for 
eadi  t.  Since  deg  F(z)  =  p  the  FIR  inverse  is  anticausal 
if  and  only  if  det  F(2)  =  cz~^  (Theorem  5.3,  [2]).  This 

will  be  the  case  if  and  only  if  UV^  has  p  eigenvalues 
equal  to  unity  (and,  of  course,  the  remaining  M  —  p 
eigenvalues  =  0).  Finally  the  FIR  inverse  is  caus^ 
(i.e.,  F(z)  unimodular)  if  and  only  if  the  determinant 
is  constant,  i.e.,  Aj  =  0  for  all  t.  W  V 

Note.  Since  UV^  has  rank  p  it  can  have  at  most  p 
nonzero  eigenvalues.  But  there  could  be  fewer,  as  in 
the  extreme  example  of  a  triangular  matrix  with  all 
diagonal  elements  equal  to  zero. 

Theorem  2.2.  The  inverse  of  the  system  F(z)  de¬ 
pends  O’)  the  properties  of  V^U  as  follows. 

1.  F(z)  has  an  anticausal  inverse  if  and  only  if  V^U 
is  nonsingular. 

2.  The  inverse  of  F(z)  is  FIR  if  and  only  if  all  eigen¬ 
values  of  V^U  are  restricted  to  be  O’s  and  I’s. 

3.  The  inverse  of  F(2)  is  FIR  and  anticausal  (i.e., 

F(z)  is  cafacafi)  if  and  only  if  V^U  has  all  eigen¬ 
values  equal  to  unity. 

4.  The  inverse  of  F(z)  is  FIR  and  causal  (i.e.,  F(z)  is 

unimodular)  if  and  only  if  has  all  eigenvalues 
equal  to  zero.  0 

Proof.  Part  1  is  a  repetition  of  Lemma  2.1.  Parts 
2  and  4  follow  from  Theorem  2.1  by  using  the  fact 
that  every  nonzero  eigenvalue  of  the  matrix  PQ  is  an 
eigenvalue  of  QP  (for  any  P,  Q  for  whidi  PQ  and 
QP  are  defined).  Part  3  follows  by  combining  parts  1 

and  2;  indeed,  nonsingularity  of  \Au  and  the  condition 
that  the  eigenvalues  be  restricted  to  be  ones  and  zeros 
is  equivalent  to  the  statement  that  all  the  eigenvalues 

of  V^U  are  equal  to  unity.  V  W 


1.  FIR  if  and  only  if  all  eigenvalues  of  UV^  are  re¬ 
stricted  to  be  O’s  and  I’s. 


Example  2.1.  The  cases  where  V^U  =  Ip  and  V^U  = 
0  give  examples  of  FIR  systems  with  anticausal  and 


2 


causal  FIR  inverses  respectively.  With  F(2)  =  (I  — 
+  z~^UV^)  we  have 


f  I  -  UV^  +  zWVt  for  Vtw  =  Ip 
1 1  +  Wvt  -  2-iWVt  for  Vti/  =  0 


as  one  can  verify  by  direct  multiplication.  Notice  that 

if  V^U  =  Ip  then  the  inverse  is  also  of  first  order.  So 
first  order  cafacafi  systems  with  higher  order  inverses 
[3]  are  not  covered  by  this  example.  Ek).  (1.1)  is  a 
special  case  of  F(z)  where  p  =  l. 

Example  2.2.  The  following  example 


UV^  = 


01 

0 

1 


0  1  0 
1  0  0 


V^Li  = 


0 

1 


0 

0 


satisfies  part  4  of  Theorem  2.2  so  that  F(z)  is  unimod- 
ular,  even  though  V^U  /  0. 

Example  2.3.  As  a  special  case  consider  I— P+2“*P 
where  P^  =  P.  With  p  denoting  the  rank  of  P,  we  can 
write  P  =  UV^  Now  P2  =  P  implies  UV^UV^  =  i/vt. 
Premultiplying  by  and  postmultiplying  with  V  and 
using  the  facts  that  U^U  and  V^V  are  nonsingular  we 

obtsdn  V^U  =  I.  FYom  part  3  of  Theorem  2.2  we  there¬ 
fore  conclude  that  there  exists  an  anticausal  FIR  in¬ 
verse  for  I  —  P  -f  2~'P,  when  P^  =  P.  In  fact  the 
inverse  is  I  -  P  -f  zP,  as  verified  by  substitution. 

Example  2.4.  Unimodular  system.  By  a  slight 
modification  of  the  above  theorem  we  can  show  that 
I  +  z~^U'\A  is  unimodular  if  and  only  if  V^U  has  all 
eigenvalues  equal  to  zero. 


3.  COMPLETE  PARAMETERIZATION  OF  BOLT 


The  BOLT  matrix  has  the  form  (1.2)  and  has  anti¬ 
causal  FIR  inverse.  Since  this  implies  G(l)  is  nonsin¬ 
gular,  we  can  write  G(z)  =  G(l)F(2l  where  F(2)  is  as 
in  (2.1).  Using  Theorem  2.2  (part  3)  we  can  say  that 
a  system  G(z)  is  BOLT  if  and  only  if  it  has  the  form 


G(2)  =  G(1)(I-WV^ -l-z-*Wvt)  (3.1) 


where  V^U  has  all  eigenvalues  equal  to  imity. 

In  Sec.  1.2  we  saw  that  the  BOLT  can  be  factorized 

as  in  (1.3)  where  V„(2)  is  as  in  (1.1)  with  vjLum  =  1 
and  p  is  the  degree  of  G(z).  Conversely  an  arbitrary 
product  of  the  form  (1.3)  with  as  above,  still 

represents  a  system  with  anticaus^  FIR  inverse,  but 
may  have  order  >  1  i.e.,  it  may  not  be  a  BOLT.  Sup¬ 
pose  we  impose  the  further  restriction  that 


l<i<fc-l, 

i  =  k 


(3.2) 


Then  it  can  be  verifed  that  (1.3)  reduces  to  (3.1)  with 
the  constant  matrices  V  and  U  given  by 

V=[vi  V2  ...  Vp],  W=[U1  U2  ...  Up], 

(3.3) 


Conversely,  can  we  represent  any  BOLT  system  as 
in  (1.3)  with  the  restriction  (3.2)?  The  answer  is  in  the 
affirmative:  if  G(z)  is  BOLT,  this  means  in  particular 
that  it  has  an  FIR  inverse,  and  so  G(l)  is  nonsingular. 
So  we  can  always  write  a  degree  p  BOLT  as  in  (3.1), 
where  U  and  V  are  M  x  p  matrices  with  rank  p.  Now 

U'iA  =  l/TT^  for  any  unitary  T,  and  we  can  rewrite 

Wvt  =  UiV\  by  defining  Wi  =  l/T  and  vj  =  Ttyt. 

Note  that  vjl/i  =  T^ytl/T.  By  proper  choice  of  T 

we  can  ensure  that  v\Ui  is  a  triangular  matrix.  In 
other  words,  we  can  assume  without  loss  of  general¬ 
ity  that  ytz/  is  triangular.  Since  G(z)  is  cafacafi  the 

matrix  has  all  eigenvalues  equal  to  unity  (part  3, 
Theorem  2.2),  so 


V^U  = 


ri 

X 

X 

...  X  ■ 

0 

1 

X 

...  X 

0 

0 

1 

...  X 

.6 

6 

6 

...  i . 

(3.4) 


where  x  stands  for  possibly  nonzero  elements.  Now 
denote  the  columns  of  V  and  ZV  as  in  (3.3).  Then  the 
property  (3.4)  means  that  (3.2)  is  satisfiea. 

Thus  we  have  defined  Um  and  v^,  1  <  m  <  p  satis¬ 
fying  (3.2).  We  already  know  that  if  such  Um  and 
are  us^  in  (1.3),  the  result  has  the  form  (3.1)  with  U 
and  V  given  by  (3.3).  So  the  given  BOLT  matrix  (3.1) 
can  indeed  be  represented  as  in  (1.3),  with  the  vectors 
satisfying  (3.2).  Summarizing,  we  have: 

Theorem  3.1.  BOLT  Characterization.  Consider 
an  M  X  M  transfer  matrix  G(z).  We  say  that  this 
is  a  BOLT  if  G(z)  =  g(0)  -I-  2“^g(l),  and  it  has  an 
anticausal  FIR  inverse.  The  following  statements  are 
equivalent: 


1.  G(z)  is  a  BOLT. 

2.  G(z)  can  be  factorized  as  in  (1.3)  where  G(l)  is 
nonsingular  and  V„,(z)  are  as  in  (1.1),  with  the 
vectors  \k  and  Ui  satisfying  (3.2). 

3.  G(z)  can  be  written  in  the  form  G(z)  =  G(l)(l  — 
IC\A  +  z“*Z/V^),  where  G(l)  is  nonsingular  and 
V^U  has  all  eigenvalues  equal  to  unity. 

4.  G(z)  can  be  written  in  the  form  G(2)  =  G(1)(I- 

+  z~^U)A),  where  G(l)  is  nonsingular  and 
\Au  has  the  triangular  form  (3.4). 


Thus  if  G(z)  is  BOLT  with  V^U  written  in  the  form 
(3.4),  the  columns  of  V  and  It  satisfy  (3.2),  and  can  be 
taken  to  be  the  vectors  v„,  and  Um  in  the  factorization 
(1.3).  So  the  factorization  is  determined  simply  by 
identifying  the  columns  of  V  and  U.  0 

Degrees  of  freedom.  The  nonsingular  matrix 
G(l)  has  elements,  and  the  vectors  Ujb,Vfc,  have 


3 


M  elements  each,  but  constrained  by  (3.2).  In  the 
real  coefficient  case  it  can  be  verified  that  the  M  x  M 
degree-p  BOLT  has  —  0.5p(p  +  1)  degrees 

of  freedom.  This  should  be  compared  with  the  LOT 
which  has  0.5M(M  —  1)  +  pAf  —  0.5p(p  +1)  degrees 
of  freedom.  IVaditional  transform  coding  (special  case 

with  =  0  and  G(l)  unitary)  has  0.5M(M  —  1) 
freedoms.  The  extra  freedom  offered  by  the  BOLT  can 
perhaps  be  exploited  to  obtain  better  attenuation  for 
the  analysis  filters,  or  to  impose  other  constraints  such 
as  linear  phase,  regularity  (for  wavelet  s^thesis  (9]) 
and  so  forth.  This  requires  detailed  investigation.  Fig. 
3.1  shows  the  magntiude  responses  of  the  8-channel 
analysis  filters  for  the  BOLT  and  the  LOT,  with  the 
same  degree  p  =  4. 


ctiam  maliix  mainx  ttiain 

Rg.  1.1  (a)  The  maximally  decimated  filter  bank,  and  (b)  the 
polyphase  representation. 


Rg.  2.1.  A  minimal  intplementation  ol  (2.1). 
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Rg.  3.1.  Magnitude  responses  ol  the  analysis  Mers 
for  (a)  the  B(XT  and  (b)  the  LOT. 
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